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Abstract 

We consider the factorisation of one-loop amplitudes at complex kinematic points. By determining 
the terms that are absent for real kinematics, we can construct a recursive ansatz for the purely 
rational pieces of one-loop amplitudes in massless theories. We illustrate this method by verifying 
the Bern et.al. n-point ansatze for the single-minus one-loop amplitudes in Yang-Mills theory and 
by constructing the scalar contribution to the one-loop five graviton MHV scattering amplitude. 
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I. INTRODUCTION 



In recent years significant progress has been made in the computation of scattering am- 
plitudes in gauge theories and gravity by utilizing the analytic properties of these ampli- 
tudes . One ingredient in this process has been the use of factorisation properties when 
the momenta in the amplitude have been complexified 1 . In particular, on-shell recursive 
methods have been very useful in the evaluation of many massless tree level processes (H, 0] • 
For example, consider shifting two of the external momenta according to: 

Aj — y Xi = Aj — z\j , . 

- _ (1-1) 
Xj — y Xj — Xj ~f" zXi 

where z is a complex parameter. Providing the shifted amplitude A(z) is analytic and 
vanishes at large \z\, then by Cauchy's theorem, we may obtain the unshifted function from 
the residues at the poles in A(z), 



A(0) 




X2) 



This is only useful if we can evaluate the residues and hence an understanding of the singular- 
ity structure of the amplitude is essential. At tree level the factorisation is relatively simple: 



amplitudes must factorise on multi-particle and collinear poles. Defining = Y^=i, kj 
as K becomes null the n-point tree amplitude A^ ee factorises as 



^tree ) ^ 



A r+ < { (/Cj, . . . , fci+r-i, K ) A n _ r+1 ( K , ki +ri . . . , 



;i.3) 



where A denotes the helicity of the intermediate state. Consequently, simple poles in the 
shifted amplitude A(z) occur at values of z where K 2 (z) = 0. Since k a + kj, is independent 
of z, only those K 7 s containing precisely one of k a or k^ will be z dependent. When the 
corresponding K 2 (z) vanishes the residue will be the product of the tree amplitudes defined at 
z = Z{. Thus the n-point tree amplitude can be expressed in terms of lower point amplitudes: 

4T(o) = E<S(^4 A ^ (Zi) ' (L4) 

where the summation over % is only over factorisations where the a and b legs are on opposite 
sides of the pole. This is the on-shell recursive expression of Q. 

There are several complications in applying these techniques beyond tree level. Firstly, 
loop amplitudes can develop higher order singularities for complex momenta . While these 
do not block recursion per se, they do necessitate an understanding of factorisation beyond 
leading order. Suppose a rational function has a double pole so that 



1 A null momentum can be represented as a pair of two component spinors — X a X a . For real momenta 
A = ±A* but for complex momenta A and A are independent 
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then 

■— + - (1-6) 

and both the leading and sub-leading, or 'pole under the double pole', terms are needed in 
order to apply recursion as in eq. f)1.2p . In general, the structure of the sub-leading poles is 
poorly understood. 

Secondly, in general loop amplitudes contain both rational and non-rational pieces. One 
strategy for computing one-loop amplitudes is to split the amplitude into a cut-constructible 
piece and a purely rational piece, 

A n = C n + R n (1.7) 
The C n may be computed using unitarity techniques 0, 0, 0- 10] and the remaining R n may 



then, in principle, be determined recursively via Cauchy's theorem provided the singularities 
of R n = A n — C n are understood. 

A general n-point one-loop amplitude in a massless theory such as gravity or QCD can be 
expanded in terms of loop momentum integrals, I m [P (£)], where m denotes the number of 
vertices in the loop and P d (£) is a polynomial of degree d in the loop momentum I. Performing 
a Passarino-Veltman reduction [TJJ on the loop momentum integrals yields an amplitude (to 
0(e) in the dimensional reduction parameter e), 

4> Ioop = E c * J 4 + E *i J 3 + E ek J * + R - > 

where q, c^-, and R n are rational functions and the J4, J3, and I2 are scalar box, triangle and 
bubble functions respectively. The mathematical form of these integral functions depends 
on whether the momenta flowing into a vertex are null (massless) or not (massive) [12]. In 
terms of this basis we can define 

c n = j2 ^1 + E ^' J 3 + E efcJ 2 (°) 

i j k 

The coefficients, Cj, dj and e^, contain a range of singularities that are not present in the 
full amplitude. Individual coefficients may contain spurious singularities of the form A~ p , 
where A is a Gram-determinant of an integral function, yet the entire amplitude is finite as 
A —7- 0. These singularities in the coefficients can be of high-order in A : we will encounter a 
case of P = 5 in one of our examples. These spurious singularities cancel amongst the terms 
in C n and also, crucially, with the rational term R n . There are also singularities that occur 
at the same kinematic points as the physical singularities, but are of higher order. Again 
cancellations between the terms in C n and R n must remove these higher order poles from the 
complete amplitude. 

Starting from C n , we can view this cancellation constraint as a means of generating parts 
of Rn as higher order poles generate rational descendants from the terms in C n . To evaluate 
the residue at a higher order pole the integral functions must be expanded to a correspond- 
ing order and the derivatives in this Taylor expansion eventually yield rational terms. In 
this way we obtain rational descendant terms whose origins lie in both the box and bubble 
integral contributions to C n . This does not completely specify R n but leaves the unspecified 
component free of these higher order singularities. 

In the following sections we describe how, by using axial gauge methods to under- 
stand the complex factorisation, we can apply recursion to the rational parts of one-loop 
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amplitudes. Specifically we consider two examples: the n-point single-minus one-loop 
Yang-Mills amplitude A n (l~,2 + ,--- ,n + ) and the five-point scalar supergravity amplitude 
M 5 (l~, 2~, 3 + , 4 + , 5 + ). The n-point amplitude A n (l~ , 2 + , • • • , n + ) vanishes at tree level and 
consequently is a purely rational one-loop amplitude. As such it has no cut constructible 
parts but it does have multiple poles in complex momentum. This amplitude was originally 
computed using off-shell methods [l3[ . In ref . [l4j a form for the sub-leading singularity was 
postulated and recursion used to (re)obtain the n-point formulae. Here we will prove, using 
axial gauge methods, the explicit form of the sub- leading term for the shift used in ref. (l4l |. 

The second example is a case where rational descendants of the integral functions in C n 
contribute to the rational terms. The example is that of one-loop five graviton scattering 
where we have a scalar circulating in the loop. This is the last remaining calculation to 
complete the five-graviton scattering amplitude. The expansion of this amplitude in the 
form of fll.Sp is plagued by high order singularities. We obtain R 5 recursively and describe 
how these high order singularities generate contributions to -R5 in addition to those coming 
directly from standard and non-standard factorisations. 



II. COMPLEX FACTORISATION 



For real momenta the factorisation of one-loop massless amplitudes is described in ref. [15 



.4 



l-loop r2 ^P 



A=± 



A r+1 p (ki, . . . , k i+r _i, K x ^j A^ r+l (( K) A , k 



u t+ri • • • 1 



ki-i) 



(2.1) 



where the one-loop 'factorisation function' F n is helicity-independent. This factorisation has 
single poles in K 2 . We refer to the singularities given in this equation as the standard factori- 
sations. Singularities not contained in eq. (12.1 j) we refer to as "non-standard" factorisations. 

For complex momenta we can acquire higher order poles. For a two-particle pole K 2 = 
(k a + kb) 2 = 2k a ■ kb = (ab) [ba\. For real momentum (a b) = ± [a b]* and so both vanish at 
the pole 2 . However for complex momenta we may have (ab) = but [ab] ^ 0. So terms 
such as [ab] 2 / (ab) 2 which are finite for real momenta can have multiple poles for complex 
momenta. These can be interpreted within eq. (12. ip as arising from the three-point one-loop 
amplitude acquiring a singularity. Specifically, the three-point all-plus (or all-minus) one-loop 
amplitude has a pole [l4T ] 



A l-loop 

A 3 



K 2 



V 1 - loop (K + ,a + ,b + ) 



(2.2) 



2 As usual we are using a spinor helicity formalism with the usual spinor products (j I) = (j \l + ) = 
u^(kj)u + (ki) and [j I] = (j + \l~) = u+(fcj)u-(fc;). In terms of spinors (ab) = e a f;X"X^ and [ab] = 
- e aB~ x tK- We also use [i\K abc \j) to denote {i + \$ abc \j + ) with.K^ bc = k%+k£+k£ etc. Also s ab = {k a +k b ) 2 , 



t abc = (k a + h + k c ) 2 , etc. 
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where, for pure Yang-Mills, 



V 1 - loop (K + ,a + ,b + ) 



487P 



[Ka] [ab] [bK] . 



(2.3) 



For real momenta A l ^ loop {K + , a + , b + ) vanishes as K 2 — > but it can be singular for com- 
plex momenta. Equation (12. 2 p specifies the double pole as K 2 — > however, as discussed 
previously, we require the subleading pole in order to use recursion. 

As an example of the structure of the double pole consider the amplitude with a single 
minus helicity A n (a~ , b + , . . . , n + ) 3 . This amplitude vanishes to all orders in perturbation 
theory in a supersymmetric theory and consequently at tree level in Yang-Mills. It is non- 
vanishing at one-loop level but, since the tree amplitude vanishes, is entirely rational. The 
all-n form was first obtained by Mahlon [13] using off-shell recursion E3- In O the complex 
factorisation of the single-minus one-loop Yang-Mills amplitude was considered by applying 
the shift of eq. (II. ip to the A of the negative helicity leg and the A of an adjacent positive 
helicity leg. For this specific case a form for the 'pole under the double pole' was proposed. 
Using this and applying complex recursion the following form for the amplitude was presented: 

^ lo °P(a-, &+,..., n + ) 

= A^(d + , . . . , n+ a" KS-^ATfr, c\ 

n-l 



be 



i=4 
i tree 



n + ,a-,K Lt )-^Ar° P (b\ 



' ' • K+J (2.4) 



+ A«™(d + ,...,n + ,a-,K b 



, V 1 ~ loop (h + r + -K + ) 

'bcJ (K^y \ U ' C ' be) 

x 1 1 + K 2 bc S^(a, K+, d)S^(c, -K~ c) b) 



where 



a, s 



b) 



(ab) 



and 



S {0 \a,s-,b) 



[ab] 



(2.5) 



(as)(sb) 1 [as][sfe] 

Expression (12. 4p was shown to match that of Mahlon [l3| up to n = 15. This expression has 
also been justified using gauge Lorentz invariance [17] . The form of the subleading pole used 
to generate ( 12 ,4p in terms of soft-factors is only valid for the particular shift used 18] . In the 



next section we provide an explicit constructive derivation of the sub-leading terms based on 
a diagrammatic analysis using axial gauge rules. 



III. ra-POINT SINGLE MINUS YANG MILLS AMPLITUDES 

In this section we study the factorisation of the single-minus Yang-Mills amplitudes 
A(a~,b + ,c + ,. . . ,n + ) under a shift of legs a and b as above. The diagrams in fig. [TJ gen- 
erate the standard factorisations given in eq. f)2.ip . As (be) — >■ double poles come from 
diagrams of the form illustrated in fig. [2] where the current r n is the sum of all possible sub- 
diagrams. To evaluate these diagrams we use axial-gauge rules Jjl. In this scheme internal 



off-shell particles are still labelled by ± helicity and the non-vanishing three-point vertices 



As usual we are considering the colour-ordered partial amplitudes. 
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FIG. 1: Standard Factorisations of A(a , 6 + , c + , . . . , n H 



are the MHV and MHV vertices 

where q is a reference null vector. For non-null momenta, P, we define 

P\q) 



\P)=P\q], \P] 



2P-q 



which corresponds to using q-nullified momenta P b , where 



P b = P 



P 2 



IP -q 



(3.1) 



(3.2) 



(3.3) 




FIG. 2: The set of diagrams corresponding to the non-standard factorisation as {be) 
gathered together as shown. r n is the set of diagrams with all but two legs on-shell. 



can be 



We can simplify this particular computation considerably by setting \ q = X a (and leaving 
X q arbitrary at this stage). With this choice three-point MHV vertices and four-point MHV 
vertices involving a~ vanish. Hence the a~ leg must be attached to a three-point MHV 
vertex and each diagram contains a single three-point MHV vertex and n — 1 three-point 
MHV vertices 4 . 



Consider constructing a one-loop diagram with n_ negative helicity external legs and n + positive helicity 
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Since the negative helicity leg a~ must be attached to the only MHV three-point vertex, 
the diagrams contributing to the non-standard (be) poles have the helicity structure shown 
in fig. [2J In this diagram r n is a current, or off-shell tree amplitude. We label the off-shell 
legs (which depend upon loop-momenta) 

B = £ + k b and C = £ - k c . (3.4) 



n 



Note that B — C = fc& + k c . As we will see, one (be) -1 factor arises from the tree current r, 
and a second from the loop integration, specifically the region where I , B and C are all close 
to null. Throughout, we view the unshifted amplitude as a sum of functions, each of which 
corresponds to a Feynman diagram involving real momenta. In particular the loop momenta 
are real and where we indicate on a diagram which legs will ultimately be shifted, that shift 
applies to the function obtained by evaluating the diagram with real momenta. 
The contribution from fig. [2] is then 

r „- sf _ f * MWj WW (Ca) 2 r n (a-, B-, -C+, , n+) 

~7 (6a) (ca) (Ba)^ £ 2 B 2 C 2 K ] 

A factor of (be)" 1 arises from the region of integration where £ 2 = 0. Specifically, since 
B 2 = £ 2 + 2£ ■ b + b 2 = £ 2 + 2£ ■ 6, around £ 2 = 0, 

We can expand r n into sub-currents which are either MHV currents or currents with a 
single minus as show in fig. El The first structure in fig. E]we label r* n . This contains an 



d+ 



VS:- 




FIG. 3: An organisation of the diagrams within r n . We have used the vanishing of the specific single 
minus current (|3.7|) . 



explicit pole and generates a further pole upon integration to give rise to the double pole 
contributions. The other structures only generate single poles and we label them r^. The 
diagrammatic expansion gives both of these contributions in terms of off-shell MHV tree 
currents: t^ hy . We can use the general results, specialised to \ q = X a , for the currents with 



external legs from three-point MHV vertices, three point MHV vertices and n.4 four-point MHV 
vertices. Then n_ = 713 + 714 and n+ = 77,3+77,4. For our situation with 77_ = 1 there is thus cither a single 
three or four point MHV vertex. 
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one off-shell momentum, P, given by [13j, [20(, 



r sm (a-,P + ,...,n + ) = (3.7) 
1+ — n+ > = - (d)(12)Tl,n)( M ) (3 ' 8) 



r MHV 



where p = ki+i + ■ • ■ + k n + fc a . Although the last expression contains an explicit P 2 factor 
in the numerator, the term in the sum with I = f contains a 1/P 2 since Pf-i = —P and so 
survives in the P 2 — > limit. We will need a simple generalisation of f l3.8|) . 



\ 2 - 2h 

r sm ( p-^ /+ ^... n +) = [ M ) r - (P - ; /+....„+) (3.10) 



This result follows diagram by diagram in r sm as only three-point MHV vertices are present 
and every one that the non-gluonic particle encounters introduces a factor of 



\ 2-2h / v 2-2/i / \ 2-2h 

(goutg) \ / [/cjPinja) \ / (PoutO) 

(PinO) J \AVovx\0) J I (Pinfl) 

the product of which gives the factor in (13.1 Op for each diagram. 
For t*" we have 




(3.11) 



n 



7^ = V 3 (P\ B\ -C-)-^r^(a-, -P~, d\..., n+) 

(Ca) 2 (a\BC\a) ^TVj (j> + C y,d + ,..., n+) 
~ (5a) 2 s bc (P, +C a> 2 



(3.12) 



This has a very simple dependence on the off-shell momenta B and C and contains an explicit 
1 / Sbc factor which, together with the pole arising from the integration, produces the double 
pole. This term also contributes to the subleading pole. 

The remaining contributions to r n arise from the second class of diagram in fig. [3J The 
integrand does not do not have an explicit pole as (be) — > and only generates a single 
pole after integration. Since this arises at C 2 = B 2 = , we can take C 2 = so that the 
r sm structures in fig. [3] have only one massive leg. In this limit we can use the formulae of 
equations (13. 8 j) and (13. 9p for currents with a single massive leg to obtain (to leading order in 
(be)) 



(Ca)' 
(Bay 



(ca)[b\l\a) (ca)(a\BKi\a) r MHV (a~, -Kr, (j + 1)+, . . . , n+) 

(cd)(de) ...(na)[ab] ~ (cd) (de) . . . (ia) (Kia) 2 K 2 



(3.13) 

where the Ki = ki+x + • • • k n + k a are fixed by momentum conservation within the t mhv 
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structures and we have made use of: 

(Ca) (ca) 
JCd) = Jafj 



+ 0((bc)) 



(3.14) 



in the relevant integration region. In this form we see that all of the contributions to C n s f 
involve the same basic integral: 



d 



, ,b\e\a)[c\£\a)[X\B\a) i (a\bc\a) 



96tt 2 (be) 



[X\2b + c\a) 



leading to 



C 



n— s f / 



-,b + ,c + ,d + ,e + ,...,n + ) 



+ 



(ca) \b\(3\a) 



n-l 

-E 



[cd)(de) . . . (ia) 



(Kia)*Kf 



where (3 = 2b + c. Setting 7 = — b so that (3 + 7 = b + c, we have 

(7a) (ba) (cd) (ca) 



(jd) (bd) (cd) (cd) 
so to leading order in (be) we have 



+ O((bo)), 



C 



n— s f / 



-,b + ,c + ,d + ,e + ,...,n + ) 



(ia)[b\(5\a) 



n-l 



(3.15) 



i [be] (a\P{b + c)\a) r MHV (a", (b + c)~, d+ ■ ■ ■ n) 
96n 2 (be) s bc (Pb+ C a) 2 ' (cd) (de) . . . (na) [a6] (3.16) 

(ca)(a|^|a) r MHV ( a -, -tfr, (i + 1)+, . . . ,n+) 



(3.17) 



i [be] (a\P{b + c)|a) r MHV (a~, (6 + c)~, rf+ • • • w) 
96tt 2 (6c) ^ (P b+C a) 2 + (ld)(de) . . . (na)[ab\ 

(a|/3^|a)( 7 a) 2 r sm (-«-, 7+, d+ • • • i+) r MHV (a', — (i + 1)+ • • • n) 



i=d 



KM 



(Kid) 



(K ia y 



(3.18) 

Again the internal momenta, «j, are specified by momentum conservation within the r sm 
factors. The quantity in the square brackets is now essentially r n (a~, (3~, 7 + , • • • ,n + ). For 
clarity we can absorb the second term into the summation by adopting an appropriate defi- 
nition for r MHV (a~, K~). 
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Using the explicit form for 

r MHV in 

the r tn contribution we have, to leading order, 

C n - si {a~,b + ,c + ,d + ,e + ,...,n + ) 



i [6c] 
96tt 2 (be) 

i=d 



(a\(3(b + c)\a) 



[qn] (a\(b + c)d\a) 



(ad)(de) ■ ■ ■ (in — l)n)(na) I [ag][an] 



tbed^bc 



n-1 

£ 

l=e 



[a\K 



a,l+l..n 



l\a) 



^a,l+l..n^a,l..n 



(a|/^|a)( 7 a> 2 r sm (- K ~, 7 +, d + r MHV (a~, -ifr, (i + 1)+ . . .„) 



As discussed previously, the contribution to the rational term is 



Res 



•s f 



V 21 



(d ,b + ,c + ,d + ,e + , . . . ,n + ) 



(3.19) 
(3.20) 



(6c) =0 



To extract the residue of the double pole term we use: 



Res 



c] 1 



1 



[ac) [b\c + d\a)(ac) 



\(bc) zt bcd s bcJ (bc)=o (bc)(a\(b + c)d\c)\(bc) (a\(b + c)d\c) 
The first term here is precisely the double pole contribution of [3] 
t(o) ,j+ „+ r >-N * T n-ioow£+ „+ r >+x_ z ( ac > (a|/3(6 + c)|a) 



(3.21) 



< u ! 1 (d+,...,n + ,a-,K, 



6c7 (^ 2 c ) 2 



r l-loop^+ )C + ; _^ 



'"■ ' " 96tt 2 (be) 2 (cd) (de) • • • ((n - l)n) (no) 

(3.22) 



The second term in f)3.2ip contains only a single factor of (be) in the denominator and its 
coefficient is unaffected by the shift. We combine this with the single pole pieces of (I3.19P to 
write the full sub-leading or pole under the pole (PUP) contribution as: 



C 



PUP 



X 



i [be] 
96vr 2 (be) 

(a\P(b + c)\a) 
(ad)(de) ■ ■ ■ ({n — l)n)(na) \ [aq][an] 



n-l 



[qn] (a\(b + c)d\a)[b\c + d\a)(ac) y-> (a\K a> i + i__ n l\ci) 



[bc](a\(b + c)d[cY 



l=e 



Sa,l-\-l..nSa,l. 



+ £ 

i=d 



(a|/^|a)( 7 a> 2 r sm (-K~, 7 +, d+ ■ ■ ■ i+) r MHV ( a -, -K~, (i + 1) + ■ ■ • n) 



KM 



(Kid)' 



(K ia y 



(3.23) 



where f denotes that the quantity in square brackets is to be shifted and evaluated at z — 
— (bc)/(ac). The sums in this expression are most of the terms in the expansion of an on- 
shell MHV amplitude. We can use the simple interchange properties of the Parke- Taylor 
amplitudes and r sm to gather many of these terms into the finite on-shell MHV amplitude 
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r MHV (a ,ft , d + , 7 + , e + • • • n + ). From a diagrammatic perspective we have 
r MHV (a-,r,rf + ,7 + ,e + ---n+) 



mhv/ — o— j+ + + +x [rf|^i|a)(/3a) 2 r MHV (a-,-^r, 7 + ,e + ---n) 



(da) (K x a) 2 K\ 

+ (a|/3/t 2 |a)(/3a)' 



v2 r sm (-« 2 -, d+, 7 + ) r MHV (a", -tf 2 " e • • • n) 



(/« 2 a) 2 «; 2 (K 2 a) 2 K 2 

+ g^ W ' ^^' , "' it)Tm( " , ^ 1)t "'" ) - (3 ' 24) 

Interchanging the 7 and d legs in the r sm amplitudes in the final sum introduces a simple 
factor of the form —((ad) / (cry)) ((7c) / (de)). We can then replace all but one term of the final 
sum in (13.231) with r MHV (a~, /3~, d + , 7 + , e + • ■ ■ n + ) and the first two terms in the expansion 
fl3T24j) : 



PUP 



X 



i [be] 
96tt 2 (be) 

(a\p(b + c)\a) 



[qn] 



(a\(b + c)d\a) [b\c + d\a) (ac) 



(ad)(de) ■ ■ ■ ({n — l)n)(na) \ [aq][an] 



[bc](a\(b + c)d\c)< 



n— 1 

E 

l=e 



[a\K 



a,l+l..n 



l\a) 



^a,l+l..n^a,l..r 



+ 



(a|/3/£ 2 |a)( 7 a) 2 r sm (- K2 -, 7 + ,d+) 



-MHV , 



a , — K 2 , e 



n 



K 2 ^ 



(n 2 a) 2 
(ad) (rye) I [d\Ki\a)(^a) 



2 MHV 



(K 2 a) 2 
'a~,-K^,-f + ,e + ■■■n) 



(aj)(de)\ (da) K 2 (K ia ) 2 

(a\(3K 2 \a)(-ia) 2 r sm (-/% , d+, 7+) r MHV (a~, -i^", e+ • • • n+) 



+ 



(jay 

(p a y 



.MHV 



(/t 2 a) 2 

(a-,r,d + ,7 + ,e + ---n + ) 



(3.25) 



The shift puts the final r MHV on-shell and we can use the Parke- Taylor form for this term. 
Using (13.91) for the off-shell r MHV factors we see that all of these contain a common term of 



[qn] 



[aq] [an] 



n-l 

E 



[a\K a . 



1 + 1..U 



l\a) 



Sa,l+l..nSa,l..n 



(3.26) 



The overall coefficient of this term vanishes before we apply the shift. The sum in (I3.9P for 
,K l ,7 + ,e + ■ ■ - n) also contains an / = 7 term. On the pole this cancels with the 
contribution from the second term in (13. 25p leaving just the r MHV (a~, d + , 7 + , e + ■ ■ ■ n + ) 
term, 



C 



PUP 



[be] 



96vr 2 (be) 



(ad)(je) (-fa) 2 
(07) (de) ((5a) 2 



_MHV 



T 



(a ,p ,d + ,7 + ,e H 



n 



(3.27) 



J t 
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Performing the shift, under which A«, A 7 — > A c at the pole, we obtain 



?pup i N „ {ad) (cc) mhv/ 



967T 2 (be) (ac) (de) 



which exactly reproduces the soft factor of ref. 14J given in eq. (12. 4p 



IV. FIVE GRAVITON AMPLITUDE 

A one-loop graviton scattering amplitude can receive contributions from a range of par- 
ticle types circulating in the loop. We denote the contribution from a particle of spin-s to 
the graviton scattering amplitude by Mn^ (with Mi ' representing a real scalar). In a super- 
gravity theory there can be contributions from minimally coupled matter multiplets. The 
contributions from the various supergravity multiplets are 2l| 5 



M^= 8 =M® + 8Mp] + 28M« + 56M^ + 70M® 

M Af=6,matter =M &/2] + g M |l] + 15A f|L/2] + 2 0M l ° ] 
M Af=4,matter =M [1] + 4M |l/2] + g M J.O] 

M Ar=l,maiter =M [l/2] + 2M [0] ( 4 _-Q 

These relations can be inverted to obtain a supersymmetric decomposition of the pure 
graviton scattering amplitude, 

jl4 2] = M^= 8 - sM^ =6 ' matter + 20M^ =4 ' ma " er - 16M^ =1 ' matter + 2M™ (4.2) 

Compared to Yang-Mills theory, the one-loop amplitudes for graviton scattering are rel- 
atively poorly understood. Previously, only for n = 4 have all the components of (14. 2 p 
been computed [21h 23(| . For n = 5, the purely rational amplitudes M 5 (+ + + + +) and 

M 5 ( h + + +) only have non- vanishing scalar components which have been computed in 

refs. [2~4j and 25| respectively. For the MHV amplitude M 5 ( h + +) only the supersym- 
metric components have been computed previously: the J\f = 8 component in ref [241 ]. the 
M = 6 in [26| and the M = 4, 1 components in [271 . 28[. In this section we use complex 



factorisation to obtain the last remaining component, M§ , of the five graviton scattering 
amplitude. 

Our starting point is (11. 8p and its the apparently trivial rewriting: 

R n = ^(Feynman Diagrams) - ^ a F 4 , trunc - ^ e k if , (4.3) 

where we have absorbed the triangle integral contributions into the truncated box contribu- 



tions [221, [29] and identified the amplitude with a sum of Feynman diagrams. The coefficients 
of the box and bubble integral functions in ( 14. 3 p are obtained using four dimensional unitar- 
ity methods and we then perform a BCFW [4J recursion on (14. 3p to obtain R 5 . Applying a 
BCFW shift 

Ai ->■ Aj = Ai - zA 3 , A 3 -»• A§ = A 3 + 2A1, (4.4) 



5 We use the normalisation for the full physical amplitudes M trcc — i(n/2) n 2 M tree , M 1Aoop 
i(27r)- 2 (K/2)"Af 1 - lo °P. 
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we obtain 



R 5 = Res(^M(^)-^Q/l :trunc -^e fc / 2 fc ) ) 

z^O poles \ ^ ieC ke£ ' / 



(4.5) 



subject to the usual caveat about the behaviour of M(z) at large z. We denote the contri- 
butions arising from the diagrams, boxes and bubbles as P diag , R hox and R huh respectively, 
so that, 

R 5 = Rf** + R h ™ + R^, (4.6) 

The rational descendants of the box and bubble contributions are obtained by expanding 
them around whatever physical or spurious singularities they contain, while the poles in 
the Feynman diagrams correspond to the standard factorisations of the amplitude and the 
non-standard factorisations discussed previously. 



A. Factorisations 



The diagrammatic contribution has two parts: the standard and non-standard factorisa- 
tions. We therefore set 

Rf a % = Rf + R™ (4.7) 



As there are non- vanishing one-loop scalar single-minus amplitudes, the five-point amplitude 
has standard factorisations (I2.ip of the form 



M tree 


[a' 


-,b~, 


-P + ) 




-loop 


M tree 


[a 


~,c + , 


~P + ] 




-loop 


M tree 


[a 


-,c + , 


-P~) 




-loop 


M tree 










-loop 



(4.* 



These contain simple poles so the rational contributions come from the rational parts of the 
four-point amplitudes. The shift (I4.4p excites six different poles, leading to 



Rf =— (M tree (l-, 2" -P+) x R 4 {P-, 3+ 4+ 5+)) 

+ — fM trce (l-, -P~, 4+) x RJ2-, P + , 3+ 5 H 
s u V 

+ — fM trce (l-, -P~, 5+) x RJ2~, P+ 3+ 4 H 
+ — fM trce (2-, 3 + , -P + ) x R 4 (i~, P-,A + , 5 H 

S23 v 

+ — (M trcc (-P-,3+,4 + ) x i? 4 (i-,2-,P + ,5 H 
+ — (M tICC (-P-, 3+ 5+) x R 4 (i~, 2", P+, 4 H 



Sl5 =0 



(4.9) 
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where the rational parts of the four-point amplitudes are 211 ] : 

„,_,+ + ,+x 1 ( s ab s ad \ 2 ( [bdf \i( 2 2 

i? 4 (a ,& + ,c + ,gT) = — ( — J ^ r ]n „ — — t ; 1 ) ( s flh + s o6 s ad + s 



[a 6] (6 c) (c d) [d a] 



i? 4 (a ,b ,c + ,d + ) 



SabSgd {a ft) 

360s^ V(6c) (cd) (da) 



' (2s 4 ad + 23s ac s 3 ad + 222s 2 ac s 2 ad + + 2s£,) 

(4.10) 



There are also non-standard factorisations. As discussed in section 3, we expect a complex 
pole when adjacent massless legs on a loop become collinear as in fig. H] 




FIG. 4: The non-standard factorisation diagram at five-point 



At five-point the tree current, r, is MHV which does not diverge in the (cd) — > limit. 
We therefore only find (cd) -1 poles in this case. The region of interest has all three of the 
illustrated propagators close to null, so to leading order we can replace r by an on-shell tree 
amplitude 

d d £ ( WM [d\£\q) \ \(e + ,A,B,a-,b-) ^ f jJ [c\e\q) [d\£\q) \ 2 M^(e+, A, B, or, b~) 



(eg) (dq) J £ 2 A 2 B 2 J V (cq) (dq) J £ 2 A 2 B 2 

(4.11) 



where A = £ + c, B = — £ + d. We use the Kawai-Lewellen-Tye relations [30( to express the 
gravity tree amplitude in terms of Yang- Mills amplitudes, 

M tree (e+, A, B, a-, b~) = s AB s ab A tiee {e + , A, B, aT, b~)A tree {e + , B, A, b~, a~) 

+ s Aa s Bb A tiee (e + , A, a', B, r)A tree (e + , a", A, b~, B) 
[eB] 4 [eA] 4 (AB) (ab) [eB] 4 [eA\ 4 (Aa) (Bb) 



(4.12) 



[eA] [Ba] [be] [eB] [BA] [Ab] [ba] [ae] [eA] [aB] [be] [ea] [aA] [Ab] [bB] [Be] 
In the region of interest (AB) ~ 0, so the first term is negligible, leading to the contribution 



ons f dH 1 [ C \£\q) 2 [d\£\q) 2 [eBf[eAf(Aa)(Bb) 

n 5-.cd J £ 2 A 2 B 2 [ ae ][be](cq) 2 (dq) 2 [aB][aA][Ab][bB] U 1 ' 

For the five-point amplitude we only need the leading term on the pole which suggests 
that the details of the off-shell continuation of £, A and B are not important. However the 
integration region may contain points where one of [aA], [bA] , [aB] or [bB] diverge. As 
the only genuine IR divergences in a diagrammatic formulation arise from propagators, we 
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rewrite all such factors in the denominator as propagators using: 

1 (Ax) (Ax) 



[xA] ' [xA](Ax) (A + x) 



+ 0((cd)), (4.14) 



where the final step involves introducing a sub- leading piece (A 2 in this case). With the 
labelling in the figure, the possible propagators involving A, B, a and b are; 

1111 

(4.15) 



(A + a) 2 ' (A + 6) 2 ' (B -a) 2 ' (B -5) 2 ' 

which fixes the ambiguity with regard to using (A ± a)~ 2 . Making use of (Ax) (By) = 
(Ay)(Bx) + 0((cd)) and setting \ q = \ a we have, 

ons f dH 1 [c\A\b) 2 [d\B\b)[d\B\a)[e\B\af[e\A\af 

5-.cd J p A 2 B 2 [ae][be](ca) 2 (da) 2 (B — a) 2 (B — b) 2 (A + a) 2 (A + b) 2 U /J 1 ' ' 

We now apply some leading order reductions: 

[c\A\b)[ba](a\A\e] = -2(6 • A)[c\Aa\e] + 2(a • A)[c\Ab\e] + 0«cd)), (4.17) 

[d\B\b)[ba](a\B\e] = -2(6 • S)[d|Sa|e] + 2(a • S)[d|S6|e] + 0«cd)), (4.18) 

leading to, 

^Tcd = C + 4°I - C - C + 0«cd» , (4.19) 

where 

d d £ 1 [c\A\b)[d\B\a)[e\B\a) 2 [e\A\a) 2 [c\Ay\e][d\Bx\e] 



/I !, 



^ 7 e 2 A 2 B 2 [ab] 2 [ae][be](ca) 2 (da) 2 (B-x) 2 (A + y) 2 

(4.20) 

The terms with x — y reduce directly to box integrals using 

(B-x) 2 (A + x) 2 = 2P~x~{(B - x) 2 ~ (A + x) 2 ) + ° i{cd)) - (421) 
Giving /IS = I® + I® + 0((cd)), with 



I [0] A 



I 



[o] 



d d £ -1 [c\A\b)[d\B\a)[e\B\a) 2 [e\A\a) 2 [c\Ax\e][d\Bx\e] 

l 2 A 2 B 2 [ab] 2 [ae}[be}(ca) 2 (da) 2 2{P ■ x) (A + x) 2 

(4.22) 

dH 1 [c\A\b)[d\B\a)[e\B\a) 2 [e\A\a) 2 [c\Ax\e}[d\Bx\e} 

x,x-.b-J £ 2 A 2 B 2 [ a 6] 2 [ae] [6e] (ca) 2 (da) 2 2(P ■ x) {B - x) 2 

(4.23) 

For x 7^ y successive reductions give 

'a=E(e,,+^U)+4 : p+o(<«i>), (4.24) 

where the terms within the sum are the box integrals given explicitly in appendix A and the 
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final term is a cubic pentagon which does not contribute to the rational term: 



AO] 

l x,x:P 



ex] 2 [ex] 2 [cx][xd][x\P\x)[x\P\x) f d d i [c\A\a)[d\B\a)(a\PxABxP\a) 



[ab] 6 (ca) 2 (da) 2 [xx] 



l 2 A 2 B 2 



(B-x) 2 {A + x) 



(4.25) 



In general, although the non-standard factorisations give both rational and transcendental 
contributions, we are only interested in the former. We take I^\-b 5 as an example: 



[ex] 2 [ex] 2 [cx}[xd}[x\P\x)[x\P\x)[x\P\a) f d d £ [c\A\a)[d\B\a)(x\BP\a) 



x,x:Bs 



[ab] 6 (ca) 2 (da) 2 [xx] 



£ 2 A 2 B 2 



(B-x) 



-. (4.26) 



We are interested in the rational term generated by the box integral of fig. [5] with the nu- 
merator given in (14.261) . 



d+ 




X 




B / 


b 








A 




c+ 


e 





FIG. 5: The box integral P^-. B 



As discussed previously, we expect this box integral to have a (cd) 1 singularity. Dropping 
the pre-factor, the I dependent part of the integral is 



rbare 

x,x:Bs 



d d l [c\A\a)[d\B\a)(x\BP\a) 



£2 A 2 B 2 



(B-x) 



(4.27) 



Appealing to (II. 8p . any divergences in this integral will be contained in the cut-constructible 
pieces so the rational terms of interest here are finite. Splitting the integration into a four 
dimensional part and an e dimensional part 31J we have, 



rbare 
x,x:B§ 



d/j, 2 J 



(// 



2 ^- l -*d A 



[c\A\a) [d\B\a)(x\BP\a) 



(£ 2 - fi 2 )(A 2 - /i 2 ) (B 2 -n 2 ) (B-x) 2 -fi 2 



(4.28) 

where all the momenta are now four dimensional. Parametrising the four dimensional loop 
momentum by 

£ = a k c + (3k d + 7 A c A d + 7A d A c (4.29) 



we see that the first three propagators of eq. (I4.28P are of the form 

Scd< 



Sd x (f(a,P, - —), 

V Srrl) 



(4.30) 
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while the final one also contains 



2B-x = as cx +(/3 -l)s dx +i[xc](dx) +~/[xd](cx) = a+j\ — )s cx + 




13-1+ 




(4.31) 



where we have absorbed a pure phase into 7'. We thus have, 



(B - xf -Li 2 = s, 



cx 



F(a,/3, 7 ',7',— , 




(4.32) 



The first three propagators of eq. (I4.28P combine with the Jacobian to give the required (cd)" 1 
pole while the fourth must introduce a spurious pole: the rational term does not involve fi, 
so we have an undetermined function depending on Sd x /s cx and s c d/s cx . As written we have 
introduced a spurious [cx] factor in the denominator, while this could be cancelled by either 
a s cx /sdx or s cx /s c d factor, either of these would introduce a different spurious factor ( [dx] 
or (cd) respectively). The rational term can therefore be no more than what is required to 
cancel the spurious poles generated by the cut-constructible pieces of the integral. 

As we are interested in the (cd) — > limit, there are two cuts we might consider: 
{d, x}{x, e, c} and {x, e}{c, d, x}. On spurious poles the coefficients of pairs of bubbles 
coincide (up to a sign) so that the logarithms in the integral functions cancel leaving a 
rational descendant. Therefore we only need to calculate one bubble coefficient, with the 
{d, x}{x, e, c} cut being the most natural. A direct parametrization yields a purely rational 
cut integral, so there are no box or triangle contributions to this cut and to leading order in 
(cd) the bubble coefficient generated by the bare integral is 



Here the [xc] -2 factor is a remnant of the spurious singularity [c\d + x\c) in the (cd) — > 
limit. Thanks to the [C|A|x) factor in the numerator the order of this spurious pole is one 
less than the loop momentum power count. 

The rational descendants of the bubbles are obtained by multiplying the bubble coefficient 
by the expansion of the difference of the two integral functions: 



h[s dx ] ~ I2M = log(— ) = M £2L±£ ^) + 0((cd)) = S —- l - S -t + \^f + ... (4.34) 

Sdx Sdx Sdx * Sj„ O Sj, 



A term involving [a;c] _1 thus produces a finite rational descendant which need not be cancelled 
by the rational term. However terms with higher order poles would contribute to a spurious 
singularity and so must be cancelled. We can read off the appropriate rational term by taking 
all the terms from the expansion that leave at least one factor of l/[scc] uncancelled. So that 
this procedure does not introduce other spurious poles we must utilize any factors of [cx] 
or [dx] present in the numerator of the bare bubble coefficient (in particular for those bare 
bubble coefficients containing [ex] factors we use [cd][ex] = [ed][cx] + [ce][cfx]). To leading 
order we can also exploit the (ad) n factor using: 




(4.33) 



s cx [cx] (xc) (da) [cx] (xd) (ca) 



+ 0((cd)) 



[cx] (ca) 



+ 0((cd)) 



(4.35) 



Sdx [dx](xd)(da) [dx](xd)(da) 



[dx](da) 
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For the B 5 contribution we note that 

= 7^> (H) \^(a\d\x] + 0([*c] -1 ) (4.36) 
and obtain the rational term, 

K h x %° B (a,b,c,d,e) =(-—(^1 i-(ad) 2 (a\d\x] I £^fR = ^"7/ R ( 1.37) 



2 

1 / [cd] \ 1 2 \ [cx](ca) 1 [cd] 2 (da) 2 (ca) 



(cd) \ [sec] / 2 I [dx](da) 2 (cd)[xc] 

The rational pieces of the other box integrals can be obtained using the same procedure, 
giving the rational term 

TZ(a, b, c, d, e) =TZ aya - A {a, b, c, d, e) + TZ b ^ A {a, b, c, d, e) + TZ a , a ; B {a, b, c, d, e) + K b)b:B {a, b, c, d, e) 

5 

+ ^2 f 7 ^^ ( a ' & ' c ' d > e ) + Kb^-.Ai (a, b, c, d, e) j 

i=l 
5 

+ ^2 ( a ' c ' e ) + ^b,a-.Bi (a, b, c, d, e) ) (4.38) 

where 'fc.j^A; denotes the full rational term generated by the box integral l x ^-.Ai (i- e - the 
bare integral multiplied by its pre-factor). There are only simple poles in these contributions 
so we find 

R? = (n(i, 2, 3,4, 5) (3 4>) I 2, 3, 5, 4) (3 6>) | ( . ^ (4.39) 
B. Box Contribution 

There is only one type of box at five-point and their coefficients are readily evaluated 
using quadruple cuts |8J: 

scalar,. - ,+ r,- ^ _ (~1) 5 (da)* (ca)V] 2 [ad] 2 [b c] (b cf (b df 

box 1 ' ' ' ! ' ]) ~ 2 (cd) 8 (be){ce}(de) ' 1 J 

Around the (cd) = pole the truncated one-mass box integral function has the expansion, 
4?runc(s, *, u, m) = u(f s log(s) + f t log(t) + f m log(m)) + u 2 f r , (4.41) 

where 

oo 

fx = ^2u j fxj(s,t,m), (4.42) 

the are rational functions and for this box U — S c di s — s ac , t = Sad and m = She- Given 
the (cd) -8 singularity in the box coefficient, on the (cd) = pole the box contributions 
produce logarithmic and rational descendants with leading singularities (cd) -7 and (cd) -6 
respectively. The logarithmic descendants combine with the bubble contributions to leave 
(cd) -1 singularities in the effective coefficients of the logarithms. The multiple poles in the 
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rational descendants cancel against terms in R 5 . 

The full rational descendant of this box contribution is 

V hox {a, b, c, d, e) = c^ lar (c + , a", d + , {b~ , e + }) x TZ 7 {s ac , s ad , s cd , s be ), (4.43) 

where, since c^ ar has poles of order eight, 1Z 7 is the expansion of u 2 f r truncated to seventh 
order in u, 

7/ . u 2 lu 3 m 1 u 4 (2s 2 + 5st + 2t 2 ) 1 u 5 m(3s 2 - 2st + 3t 2 ) 

H 7 (s, t, u, m) = -^r — H — (4.44) 

v ' ' ' ; s 2 t 2 3 sH 3 12 sH 4 30 s 5 t 5 v ; 

1 u 6 (48s 4 + 84s 3 t-40s 2 t 2 + 84st 3 + 48t 4 ) 1 u 7 m(10s 4 - 8sH + 9s 2 t 2 - 8st 3 + Wt 4 ) 



720 sH 6 210 s 7 t 7 

The shift excites the (cd) = type pole in four of the six box contributions, leading to 

R b 5 ox = Res ( 3 4)=0 l -fp box (i,2,3,4,5) + P box (2,l,3,4,5)) I (1.1 r ) I 



z 



+ Res ( 3 5)=0 ( - z (V h0x (l, 2, 3, 5, 4) + V hox (2, 1, 3, 5, 4)^ 



Bubble Contribution 



The bubble coefficients can be obtained from two particle cuts using canonical forms [10 . 
There is a single type of bubble with coefficient c b ^ ar ({a _ , c + }, {b~ , d + , e + }). The general 
form of the bubble coefficients is given in appendix [Bj 

The bubble coefficients contain poles of the form (cd)~ 7 . On the (cd) — > poles these terms 
are precisely those required to cancel the multiple poles in the coefficients of the logarithmic 
descendants of the boxes. The bubble coefficients also contain spurious poles of the form 
[d\a + c\d)~ 5 (two powers worse than the M = 1 case). On these poles pairs of bubble 
contributions combine to produce rational descendants with [<i|a + c|<i)~ 4 spurious poles. The 
singular pieces of these rational descendants are 

P bub (l, 2, 3, 4, 5) = (<S(1, 2, 3, 4, 5) + P {3 ,4,5}) + (5(2, 1, 3, 4, 5) + P {3 ,4,5}) (4-46) 
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where ^{3,4,5} represents permutations of the positive helicity legs and 



S(a, b, c, d, e) 



(ac) 4 (bd) 6 [b e] s 



cd 



60 (be) (cdf (dey[d\a + c\d) 



x 



[d|a + c|d) 3 2 [d\a + c\d} 2 [s 



s,d + 3f 1 + 3 ( cb )( ad ) (ce)(ad) 



+ 



Scd 



12 [d\a + c\d) 



1 + 3- 



ac) (6 c?) (ca) (de) 
cb) (ad) (ce) (ad) 



+ 30 3 



(cb) 2 (ad) 2 
(a c) 2 (b d) 2 



[ac) (bd) (ca) (de) 

(cb) (ad) (ce) (ad) (ce) 2 (ad) 2 
(ac) (bd) (ca) (de) + (ca) 2 (de) 2 



+ 



2*L 



+ 5 



^ _ 3 (cb) (oj) \ _ 15 s c rf (ac?) 2 (cb) (cd) (be) ^ (ce) (ac?) 
(ac) (bd)J 4 s ac (ac) 2 (bd) 2 (de) (ac) (de) 



(c d) 2 f 3 (abf_ + 1 (ae) 2 + g (a 6) (a e) 



(ac) 2 V (bo 1 ) 2 (de) 2 (bd) (de) J (ac) 3 (bd) 2 \(bd) (de) J _ 



(c d) 3 (a b) 2 f (ab) (a e) 
+ 3- 



(4.47) 



The shift excites four distinct spurious poles: [3|1 + 4|3) = 0, [3|1 + 5|3) = 0, [4|1 + 5|4) = 
and [5|1 + 4|5) = and we have 



1,2,3,4,5)= Res 



[3|i+4|3>=0 



-P bub (l,2,3,4,5) +Res [ 3 |i 



[3|1+5|3)=0 



-P bub (l,2,3,4,5) 



+ R es [4|i+5|4>=o I "^bubll, 2, 3, 4, 5) \ + ReS[ 5 |j +4 | 5 ^ =0 -P bub (l,2,3,4,5) 

(4.48) 

The residues of all terms are readily extracted using Mathematica. This completes the 
computation of -R5. Technically, the validity of the above result relies upon, a priori, the 
vanishing of R(z) for large z for which no general theorems are available. Our expression 
satisfies several non-trivial consistency conditions which, experience suggests, make it al- 
most certainly correct. In particular, the resulting expression for the amplitude has the 
correct symmetries, is free from spurious poles, has the correct soft and collinear limits and 
has the correct complex factorisation. The Mathematica expression for R$ is available at 



http : //pyweb . swan. ac . uk/~dunbar/gr avit on . html 



Our results are consistent with the suggestion that gravity perturbation theor y h as a 
significantly softer ultra-violet behaviour at one-loop than traditionally expected [26, 32[ 33| . 
If a Feynman diagram has n-points in the loop and has a loop momentum polynomial of 
degree m then bubble coefficients have (a • P) spurious singularities of the form 



1 



(a-P) 



m— n+3 



(4.49) 



and the rational terms have one power less. For gravity the traditional expectation of m = 2n 
leads to (a-P)~ n ~ 2 . However, our explicit calculations indicate a softer behaviour. The bubble 
coefficients (IB11I) have explicit (a- P)~ 6 singularities arising from the H\ x terms. For the five- 
point case the leading singularity vanishes leaving (a ■ P) -5 singularities - and consequently 
only (a ■ P) -4 singularities in the rational term. These singularities are thus consistent with 
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an effective power count of m = n + 4. This is in agreement with the expectations of 26|, [33 . 



V. CONCLUSIONS 

Recursive techniques based on the factorisation of rational tree amplitudes have proved 
of great utility. In extending these methods to one-loop amplitudes one faces obstacles 
of various types: in general one-loop amplitudes contain both rational and transcendental 
functions along with a plethora of spurious poles and higher order physical poles. The 
transcendental pieces of one-loop amplitudes are readily obtained using four dimensional 
unitarity techniques. When the coefficients of these transcendental functions contain high 
order poles these contributions give rational descendants which must be accounted for when 
computing the remaining rational terms recursively. Further, using recursive techniques for 
the computation of one-loop amplitudes requires an understanding of the singularities of 
these amplitudes when extended to complex momenta. For complex momenta there are 
"non-standard" factorisations which must be accounted for. 

In this article we have demonstrated how axial gauge techniques may be used to de- 
termine the non-standard factorisations. Axial gauge techniques provide a natural method 
for examining complex structures since they preserve the language and structure of on-shell 
amplitudes, however, like usual Feynman diagram techniques, they can prove rather cumber- 
some. Nonetheless, we have used these techniques to A) prove the conjectured sub-leading 
pole in the single-minus Yang-Mills amplitude, B) compute the rational terms of the five 
graviton MHV amplitude analytically. This completes the calculation of the five graviton 
scattering amplitude. The expression we obtain is not particularly simple but will provide a 
significant target for alternative techniques. Our results are consistent with the suggestion 
that gravity perturbation theory has a significantly softer UV behaviour than traditionally 
expected. 

The work of S. Alston was supported by a STFC studentship. 



Appendix A: Box Integrals Contributing To The Non-standard Factorisations 

The non-standard factorisation term can be expressed as a sum of box integrals: 

4S = + zEU) + i%. P + o((cd}) . (ai) 
i=i 

Using the ^-dependent variables A = £ + c, B = —£ + d, the integrals required are 
[ex] f d d £ [c\A\a)[d\B\a)[e\B\a)[e\A\a)[c\A\x)[d\B\x)[e\A\x) 



S:Al [ab] 2 (ca) 2 (da) 2 [xx] J £ 2 A 2 B 2 (A + x) 2 

[ex] f d d £ [c\A\a)[d\B\a)[e\B\a)[e\A\a)[c\A\x)[d\B\x)[e\B\x) 



1 [ab] 2 {ca) 2 (da) 2 [xx] J £ 2 A 2 B 2 (B - x) 



(A2) 



(A3) 



[o] [ex] [xe] [cx] f d d £ [c\A\a)[d\B\a)[e\B\a)[e\A\a){x\AB\x)[d\A\x) 

X ^ M ~ [ab] A (ca) 2 (da) 2 J £ 2 A 2 B 2 (A + x) 2 ' ^ ' 
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loj [ez] [are] [car] F d d £ [c\A\a)[d\B\a)[e\B\a)[e\A\a)(x\AB\x)[d\B\x) 

X '*- B2 ~ \abY(ca) 2 (da) 2 J £ 2 A 2 B 2 (B - x) 2 ' ' ' ^ 



[ex] 2 [xe][cx][xrf] F d d £ [c\A\x)[d\B\x}[e\A\x}(a\AB\a) 2 

'x.x:A : < ~ "..7.1 i /, ..\-">/./..VH- ..i / /■■> r> ;>■'»' i T~^\2 ' l AD / 

\ l„\ A Dl„\2 

(A7) 



^ :A3 [a&] 4 (ca) 2 (da) 2 [xx] J £ 2 A 2 P 2 (A + x) 

[o] [ex] 2 [xe][cx][xd] y [c|A|x)[d|P|x)[e|P|x)(a|AB|a) 2 



X ' S]B3 [ab} 4 (ca) 2 (da) 2 [xx} J i 2 A 2 B 2 [B — x) 2 

[o] [ex] 2 [ex] 2 [cx][xrf][x|P|x) /" d d £ [c\A\a)[d\B\a)(x\AP\x)(a\AB\a) 

x ' x:Ai ~ [ab] 6 (ca) 2 (da) 2 J i 2 A 2 B 2 (A + x) 2 ' ( 



[u [ex] 2 [ex] 2 [cx][xd][x\P\x) F d d i [c\A\a)[d\B\a)(x\BP\x)(a\AB\a) 

x ' x:B * ~ [ab] e (ca) 2 (da) 2 J l 2 A 2 B 2 (B - xf ' (: 9j 



,,„ [ex] 2 [ex] 2 [cx][xd][x\P\x)[x\P\x)[x\P\a) F d d t [c\A\a)[d\B\a) (x\AP\a) 



x ' x:A5 ' [abf(ca) 2 (da) 2 [xx] J £ 2 A 2 B 2 (A + x 

[o] [ex] 2 [ex] 2 [cx][xc?][x|P|a;)[a;|P|x)[a;|P|a) F d d i [c\A\a)[d\B\a){x\BP\a) 

X ' X]B5 = [abf(ca) 2 (da) 2 [xx\ J t 2 A 2 B 2 [B - x) 2 



(All) 



Appendix B: Bubble Coefficients for n-point MHV amplitudes 



In this appendix we present a general formula for the bubble coefficients of the n-graviton 
MHV scattering amplitude. The bubble coefficients vanish for M = 8, 6 contributions. The 



M = 4 contributions were given in ref. 26j and the Af = 1 coefficients in ref. 27|. We 



present the remaining scalar contribution together with these in a unified way using canonical 
forms 1(|. In the method of canonical forms the bubble coefficients are obtained from 
unitarity by decomposing the product of tree amplitudes appearing in a two-particle cut into 
canonical forms 



^M trcc H 



x M tree f- 



/* ] c iJ~i {£j 



(Bl) 



where the c, are coefficients independent of £j. The momentum across the cut is P = i\ — £2- 
We then use substitution rules to replace the J-%{£j) by the canonical forms Pj(P) to obtain 
the coefficient of the bubble integral ^(P 2 ) as 



(B2) 
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The general types of canonical forms we need for MHV amplitudes are 

n 

n^llwe^Bi) (B3) 



i=i 

n 



= JIIAI^IBi)^^. (B4) 

n^flm.m^^ <B5) 



1=1 



1=1 



Supersymmetric amplitudes require canonical forms with lower values of n. In Yang-Mills 
theory we need values n < 2 in general but only n = for supersymmetric contributions. For 
gravity, scalar contributions require canonical forms up to n = 4 with only n < 2 for Af = 1 
and n = for Af = 4. The expressions for the canonical forms with n > 2 were not given 
previously. 

We first note some general results (not given previously) 



1 n 
#o {A}] = — -ry £ II[^I P I^> 



(n+ 1)! 

i " (p2\n-r r n ~ r 

^;W;{A}] = 7 ? — E E U P \A)n-r + i UW P \ *) UlDilAlBMAlPlBn 

VV ''' P{Bi),P{Di) r=0 L I I / j =1 J=1 

(B7) 
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H\ x and H\ x are given in ref [10|. The expressions for H\ x and H^x are 
Hl[A; {B u B 2 , B 3 , B,}; C l5 {£> x , D 2 , D 3 }; P] = £ 

P({fl i }),P({£» i }) 



(P 



2\3 



[A\P\AY 



4[ J D 1 |A| J B 1 )[ J D 2 |A| J B 2 )[ J D3|^|53)[^|P|54)[A| J P|C 1 ) 

-4[ J D 1 |A|fi 1 )[ J D 2 |A|fi 2 )[ J D 3 |A|C 1 )[A|P|i?3)^|P|i?4) 



+ I 2[D 1 |A|B 1 >[D 2 |P|B 2 >[D 3 |P|C 1 >[A|P|B 3 >[^|i'|54> 



[A|P|A>* 



4[ J D 1 |A|C 1 }[ J D 2 |P|Pi)[P3|P|S2)[^|P|S3)[^|P|P4) 



+ [A| p |A)2 f [£>i|P|Si)[P> 2 |P|S 2 )[£>3|P|C' 1 )[A|P|S3>[A|P|S4)j ] (B* 
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and 



H* X [A; {B u P 2 , P 3 , P 4 , P 5 }; C i; {D u D 2 , P 3 , P 4 }; P\ = — £ 

P({B i }),P({£» i }) 

-2[ J D 1 |A|fi 1 )[ J D 2 |A|5 2 )[ J D 3 |A|fi3}[ J D4|A|C 1 )[A|P| J B 4 )[A|P|fi 5 ) 



H|P|A)' 



-22[P 1 |A|P 1 )[P 2 |A|P 2 )[P 3 |^|C 1 )[P4| J P| J B 3 )[A|P|P 4 )[A|P|P 5 ) 



+ ri.CiL 1 SS^lAlP^^l^lP^^lPiP^^lPlP^tAiPiP^tAlPlCO ] 



[A|P|A)< 
[A|P|^ 



+ r JpiL l 27[D 1 |A|B 1 )[D 2 |P|B 2 )[£)3|P|B3)[^4|i J |Ci)[A|P|B4)[A|P|B 5 ) 



-45[ J D 1 |A|C 1 )[P> 2 |P|P 1 )[P>3|P|P 2 )[P> 4 |P|P 3 )[^|^|S4)[A|P|P 5 ) 



+ TTHl^l -15[^i|P|Pi)[^ 2 |P|P 2 }[^3| J P|53)[^4|P|P4)[^|P|5 5 )[A|P|C 1 ) 



[A|P|A>* 



+ 



21[D 1 |P|P 1 )[D 2 |P|P 2 }[D3| J P|53}[^4|P|C 1 )[A|P|P4)[^|P|5 5 ) j 



(B9) 



We do not need to present P{\ separately since 

Hi tl [Ai, A 2 ; Pi, • • -P n +i; Ci; Pi, • • • D n ; P] = H 2 [Ai, A 2 ; Pi, • ■ • P n +i, Ci; Pi, • • • P] 
+ (CiA 2 )(A 2 Pi) g „_ 1[ ^ ; ^ ^ ^ P|Pi]; P 2 , • • • P n ; P] 

_ (CiA 2 )(PiAi) 1[Ai , Aa . ^ ^ p|£)i] . ^ A . p] 

(BIO) 

We can now use these canonical forms to provide expressions for the coefficients of the 
bubble integrals I 2 (P 2 ). The bubble integral functions h(P 2 ) will have vanishing coefficients 
for the MHV amplitude unless the momentum P (and hence — P) contains exactly one 
negative helicity leg and at least one positive helicity leg. We can thus take P to be of the form 
{m^, af, a 2 , • • • , a+ L } and the legs on the other side of the cut to be {m^, &+, b 2 , ■ • • , b+ R }- 
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Using the canonical forms we find the bubble coefficient 

(mi m 2 ) M 



c(mi, R}; m 2 , {fy}) - , p2 , 4 _ M 



xe{ai}U{fej}— ai 

(ai 61) 



2j Cp^CpU 

{h ai) 



E A^^Af^h; *; {A}; mi; {A}; P] 



are{ai}U{6j}-6i 



+D ai ^^-H^ M [a i; {A}; m i; {A}; P] + ^^^ M [&i; {A}; m i; {A}; P] 

(Bll) 

where Pj, and Pp are permutations of the positive helicity legs {a^} and {bi} respectively. 
The different cases are specified by M = 2 for the M = 4 matter multiplet, M = 1 for the 
M = 1 matter multiplet and M = for the scalar contribution. The arguments are 

{Bi] N=4: = {mi}, {Bi} N=1 = {mi, ma, m x }, {A} [0] = {m x , m 2 , m 1 ,m 2 , mj 
{A}^ =4 = 0,{A} A ' =1 = {P|mi),P|m 2 )},{A} [0] = {P|mi),P|m 2 ),P|mi),P|m 2 )} (B12) 

and 

C Pl = n — 1 ' ^Pfl = n — — 1 ' (B13) 

(n L mi)nr=i ( a i a i+i) " ( n R m 2)Hj=i ^+i> 
(mix) u:^ l HKi + i\x) n^r'^i^+ik) _ mr'kixi+ik) n££iM^ + ii*> 



An 



(B14) 

where K p = k ap + . . . + k mi and A^ = kb p + ■ ■ ■ k bnR + k m2 . The bubble coefficients satisfy, 
nontrivially, the IR relation ^\ q = 34J . 
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